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TECHNICAL NOTE 3451 


ANALYSIS OF lULLY DEVELOPED TURBULENT HEAT TRANSFER AND FLOW 
IN AN ANNULUS WITH VARIOUS ECCENTRICITIES 
By RoBert G. Deissler and Maynard F. Taylor 


SUMMARY 

A previous analysis for t\arbulent heat transfer and flow in tubes 
was generalized and applied to an a.nrm lns with various eccentricities. 
Expressions for eddy diffusivity which were verified for flow and heat 
transfer in tubes were assmed to apply in general along lines normal 
to a wa3J-. Velocity distributions, wall shear-stress distributions, 
and fhictlon factors, as well as wall heat-transfer distributions, wall 
tenperature distributions, and average heat-transfer coefficients, were 
calculated for an annuliis having a radius ratio of 3.5 at various 
eccentricities . 


INTRODUCTION 

Most of the existing analyses for txirbulent flow and heat transfer 
in passages have been confined to circtilar tubes or parallel plates 
(refs. 1 to 5). Th.ese passages have been analyzed extensively because 
of their importance in technical applications and because their simplic- 
ity makes them amenable to analysis. 

In recent years the problems associated with the use of odd- shaped 
passages in heat exchangers have become important. In reference 6, tem- 
perature distributions for rectangular and triangular ducts were calcu- 
lated by using experimental velocity distributions and average heat- 
transfer coefficients. No attempt was made to calculate the velocity 
distributions or heat-transfer coefficients. Seme work on the calcula- 
tion of the velocity and shear-stress distributions in comers is re- 
ported in reference 7. 

As a part of a general investigation of heat transfer and flow in 
passages of various shapes being conducted at the NACA Lewis laboratory, 
the eccentric annulus was analyzed. The methods used in references 4 
and 5 for calculating the heat transfer and friction in tubes were ex- 
tended and applied to annuli having various eccentricities. The methods 
used herein should be general enou^ to be applied to passages having 
various other cross sections. 
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BASIC EQUATIONS AND ASSUMPTIONS 

The differential equations for shear stress and heat transfer can 
be vTritten in the following form: 


T 




+ p,) ^ 


( 1 ) 


q = - (k + pgCp€h) ^ 


( 2 ) 


where e and are the eddy diffusivities for momentum and heat 

transfer, respectively, the values for which are dependent t;ipon the 
amount and kind of turbulent mixing at a point. In these equations 
y is taken as the peipjendiculax distance from the wall. Equations 
(l) and (2) can be considered as definitions of e and They 

can be written in dimensionless form as 


-I- 3 /iL + _£_ f <3u+ 

To“^4o Po'i^o7^)dy+ 

= /JL -L. + -E_ a ^ ^ 

qo \^kQ Pro PQ Cp^ Hq/pqJ dy+ 


where the subscript 0 refers to values at a wall. All symbols are de- 
fined in appendix B. 


Expressions for eddy diffusivity . - In order to use equations (3) 
and (4), the eddy diffusivity € must be eval^uated for each portion of 
the flow. It is assumed, as in references 4 and 5, that in the region 
at a distance from the wall the mechanism for turbulent transfer is de- 
pendent only on the velocities in the vicinity of the point measured 
relative to the velocity at the point or on the space derivatives of 
the velocity. In the region close to the wall the turbulence is assumed 
dependent on qziantities measured relative to the wall, that is, u and 
y.l Dimensional analysis then gives, for the region close to the wall 
(y+ < 26), 


e = n2uy 


(5) 


where the constant n has the experimentally determined value 0.109 
(ref. 4). 


lit was found in reference 8 that in the region very close to the 
wall e appears to be a function also of kinematic viscosity, but the 
effect of that factor becomes important only at Prandtl numbers appre- 
ciably greater than 1.0. 
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In the region at a distance from the -wall (y'*' > 26), e is assumed 
to he dependent on the relative velocities in the neighborhood of the 
point. From a Taylor's series expansion for u as a function of y 
and z, 

^ ~ Sy2" ■ ■ ■ ^zSy 

where y and z are measured in normal directions in the cross sec- 
tion of the passage. In the case of flow in a tube or between par’allel 
plates, the derivatives in the z-dlrection are zero because the velocity- 
gradient lines are straight lines normal to the wall. (A velocity- 
gradient line is a line which at each point is normal to a constant - 
velocity line.) In an eccentric annulus the velocity-gradient 3j.nes 
near the wall also are normal to the surface, but they are usually 
curved in the center portion of the passage. Inasmuch as the greatest 
changes of velocity with respect to distance take place in layers near 
the wall, the effect of the derivatives with respect to z will be neg- 
lected. It seems reasonable to expect that near the center of the flow 
passage the effect of the derivatives with respect to z would be to 
increase the turbialence and flatten the profile in that region. However, 
the normal t'urbulent profile (derivatives with respect to z absent) is 
already very flat in that region, so that the increased t\arbulence should 
not produce significant changes in the values of the velocities. There- 
fore, the expression for e for y^ > 26, obtained by using dimensional 
analysis, and with only the first two derivatives with respect to y con- 
sidered, is 

(d^u/dy^)2 

where x has the experimental value 0.36 (ref. 4). Hils is von Easmian's 
expression (ref. l) . 

Further assunptlons . - In order to integrate equations (s) and (4), 
the following assumptions are made in addition to those concerning the 
eddy dlffusivlty (eqe. (5) and (6)): 

(1) The fluid properties can be considered constant and the Prandtl 
number close to 1.0 (0.73). The heat-transfer results are therefore ap- 
plicable to gases with moderate temperature differences. The analysis 
could be carried out for variable properties, but the complexity would 
be increased. 

( 2 ) The eddy dlffusivitleB for momentum e and heat transfer 

are equal, or a = 1. Previous analyses for flow in tubes based on this 
assumption yielded heat-transfer coefficients that agree with experiment 
(refs. 2 and 5). At low Reynolds or Peclet numbers (Pe = Re Pr), a 
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appears to "be a function of Peclet number (ref. 9), but for Reynolds 
nimibers above 15^000, as in the present analysis ^ a is nearly constant 
for gases. 

( 3 ) Along the lines normal to the -wall, the variations of the shear 
stress T and heat transfer per unit area q have a negligible effect 
on the velocity and temperatirre distributions. It is shown in figure 11 
of reference 5 that the assvanption of a linear variation of shear stress 
and heat transfer across a tube (t or q = 0 at tube center) gives 
very nearly the same velocity and temperatiire profiles as those obtained 
by assuming uniform shear stress and heat transfer across the tube . 

( 4 ) The molecular shear-stress and heat-transfer terms in equations 
( 3 ) and ( 4 ) can be neglected in the region at a distance from the wall 
(y+ > 26)(ref. 5, fig. 12). 

Generalized velocity and teaperattare distributions . - Equations ( 3 ) 
and ( 4 ) have been integrated using the foregoing assumptions in refer- 
ences 4 and 5, Tdiere the equations are also verified experimentally for 
flow and heat transfer in tubes. The reB^Ilts are reproduced in figures 
1 and 2. These curves give the relations between u+, y+^ and t+ 
which will be used in the following calculations for eccentric annuli. 


GALCULATIOT OP -VELOCITr DISTRIBUTIONS FOR ECCEHTRIO ANNULI 

For applying the relation between u+ and y+ in figure 1 to the 
cadculatlon of velocity distributions in an eccentric annulus, an iter- 
ative procedure must be used, inasmuch as the lines of velocity gradient 
(lines normal to constant-velocity lines) are unknown at the outset. An 
eccentric einnulus, with several assimied velocity-gradient lines, is 
shown in sketch (a) . 
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By using these assumed velocity-gradient lines, lines of constant veloc- 
ity can he calculated, as will be shown. A new and more accurate set 
of velocity-gradient lines can then be drawn, inasmuch as they must be 
normal to the constant -velocity lines . With a little practice, the 
velocity-gradient lines can be estimated quite acc^r^ately the first 
time so that it is not usually necessary to go through the procedarres 
more than once or twice. 

Calculation of line of maximum velmcities . - After the velocity- 
gradient lines have been estimated, iiie next step is to calculate the 
line of maximum velocities shown on sketch (a) . The velocities on ei- 
ther side of this line are lower than the velocities on the line. Two 
adjacent velocity-gradient lines are also shown in the sketch. The 
lengths A2i and AZ 2 are on the inner and outer cylinders, respec- 
tively, and the line dividing AAq and AA 2 is the line of maximum 
velocities . The locations of AI^ a.nd AI 2 must be' chosen so that 
straight lines drawn normal to them at their midpoints (shown dashed 
in the sketch) meet at the line of maximum velocities . It will be 
necessary to match the velocities at -that point as calculated tram. 
the relation between u+ and y+ along the normals from AZq anil 
AZ 2 - In order to 'write force balances on the elements AA^ nnfl M 2 , 
the shear forces acting on AZj_ and AZ 2 and the pressure forces 
acting on the faces of the elements must be considered. There are no 
shear forces acting on the velocity-gradient lines, because the normal 
velocity deri'vatives are zero along those lines. Writing force bal- 
ances on Mj^ and M 2 results in 

'1 “ ■ ^ (®) 

^ ( 8 ) <•) 

where the pressure gradient dp/dx is uniform over the annulus because 
the flow is fully developed. Dividing equation (8) by equation (?) 
gives 


T2 _ ^2 

Ti ^^2 "^1 


From the curve of n't against y+ (fig. l). 


( 9 ) 


k 


u+ = F(y+) 
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Applying this equation to AAj_ and AA.£ gives 


n-, 


l/p 



( 10 ) 


and 






= F 



( 11 ) 


Wtien these equations are evaluated at = Ug = Uj^^, and equation (lO) 
is divided by equation (ll ) , there results 




l/P , 

ix/p ^2 m} 


%/p 


( 12 ) 


F 


Elimination of and Tg from equation (l2) by use of equations (?) 

to (9) and conversion to dimensionless form give 


■where 


^1 ^ 



_ 

^ 


u+ 


is a type of Reynolds mam- 


ber -which is assigned an arbitrary value. The quantity y-, m/rj ^ > that 

is f a point on the line of maximum velocities , can be calculated from 
equation (13) by trial. Referring to sketch (a), the -value for 

is first assumed. Values for , and y^^/r-i 

can then be calculated by measuring the areas and lengths in the sketch. 


yin 
^iuM 


and y ^ (terms in parentheses, 

■gjjj obtained from the plot of the 


For a given -value of rj+, -fchen, 
eq. (13)) can be calciilated and 
relation between u+ and y+ (fig. l) . If this value of u|^/u^ does 

not agree with the term on the extreme left side of equation (l3), an- 
other -value for y-,^/r-, must be tried. In this way -the value for 

^ thus a point on the line of maximum velocities, can be found. 
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Calculation of lines of constant velocity . - For calculating veloc- 
ities at various points In the annulus, it Is convenient to define the 
follo'wing velocity p^ameter: 



Ui 




dp/dx 


up- = 


Uc 




rj_ dp/dx 


(14) 


This pajfameter Is used In place of u"*", because the shear stress In the 
definition of u+ varies with position. Equation (l4) can be written 
In terms of quantities already known as 


where 



(15) 


(16) 


where the function F Is obtained from the relation between u"*" and 
y+ in figure 1. Equations (l5) and (l6) apply to points lying between 
the inner cylinder and the line of maxlnnjm velocities. For points be- 
tween the outer cylinder and the line of maximum velocities, these- equa- 
tions are replaced by 


where 



(17) 


(18) 


The quantities y^, y-,^/r.j , and y ^ /r^ are already known from 

equation (l3). The relation between u^ and y-^/y.,^ or u^ and 
yg/ygm s- given r^ can therefore be calculated along a given 

line normal to the inner or outer cylinder. 
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By carrying out the calculation for various lines normal to the 
inner or outer cylinder, lines of constant velocity (constant u++) cam 
be obtained. As mentioned previously, new and more accurate velocity- 
gradient lines are next drawn so as to intersect the constant-velocity 
lines at right angles. The calculation is then repeated using the new 
velocity-gradient lines. 

Calculated velocity distributions . - Several velocity distributions 
for eccentric annuli, as calciolated by the method described, are shown 
in figure 3. The calculations were carried out for an annulus having a 

diameter ratio of 3^ for various eccentricities. The shape of the 

constant-velocity lines indicates, as would be expected, that the great- 
est velocities occur on the side 'vdiere the separation of the cylinders 
is greatest, and that the velocities on the side where the separation is 
least are much smaller and go to zero idien the cylinders touch. Although 
the velocity distributions are plotted for an rj^ of 200, the constant- 

velocity lines for an r^ of 4000 have practically the same shape. The 
values of are slightly different for the two values of r^, 

but the difference is not large. 

In 1 cases the line of maximum velocities lies closer to the 
inner than to the outer cylinder, that is, it lies closer to the stir- 
face having the smaller area. 

In most cases the calculated constant -velocity lines shown in fig- 
ure 3 are very nearly normal to the velocity-gradient lines, as required. 
However, in some cases, especially for the large eccentricities, diffi- 
culty was experienced in obtaining noimal lines in some regions. Whether 
this lack of normalcy is due to the fact that the iterations were not 
carried far enough or to the approx±raate natinre of some of the assump- 
tions made in the analysis has not been established. However, the diffi- 
ciilty occurs only in regions in which the velocity gradients are very 
small, so that the error in the computed values of velocity should be 
small. 


lEICTIOW I’ACTOR AND EEYHOLDS NUMBER 

Once the velocity distributions for the annuli have been obtained, 
friction factors and Reynolds numbers can be calculated by integrating 
the distributions to obtain bulk or average velocities. The bulk ve- 
locity between two adjacent straight lines normal to the wall (dashed 
lines in sketch (a)) is first obtained. This bulk velocity varies with 
position on the inner cylinder and is given by 
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(19) 


■where AA is the to’tal area between t'wo adjacent straight lines normal 
to •the wall (on both sides of the line of Tnax1.ntuTn velocities). Equation 
(l9) can be -written in dimensionless form as 



( 20 ) 


The average bulk velocity for the whole annulus can be written as 



( 21 ) 


where A 

“T/“btav 


is the total area of the annulus. The variation of 



■with angle and eccentricity is shown in figure 4. 


For later comparison with heat transfer, the local friction factor 
based on the shear stress on the inner cylinder and the average bulk ve- 
locity is introduced and is defined as 


f = 

P^b,av^ 

or, in tenns of known dimensionless quantities, can be -written as 

r - / T 

■^1 l^btav ^l*'(yim/^l)J 

The average friction factor based on the shear stress is, then. 


( 22 ) 


(23) 


^l,av 


1 

3t 



f^ a.0 


(24) 


where 9 is the angle at the center of the inner cylinder measured from 
the line of least separation of the two cylinders. Di-vision of equation 
( 22 ) by equation (24) gives 
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fT-. 




”^l,av 


(25) 


The ratio have heen obtained directly from eqmtion (?) 

and a corresponding equation for the whole annulus as 


"^l^av ” -^1 -^1 

Curves for Tt/tt „„ or f, /f_ against angle 0 are pre- 

j. x,ttv ‘i^av 

seated in fig^nre 5'. The cvirves indicate a decrease in shear stress on 
the side of the inner cylinder \dalch is closest to the outer cylinder, 
and a shear stress of zero where the cylinders touch (e = l.O). As in 
the case of the velocity distributions (fig. 4) the value of rj^ or 

Reynolds number has a negligible effect on the shape of the curves. The 
values for the Reynolds ntmbers were calculated from 


Re 


Ui 


++ 


^av 


We 




(26) 


which follows directly from the definition of Reynolds number. The 
eq,uivalent diameter Dg is equal to the difference between the diame- 
ters of the outer and inner cylinders. A cross plot of figure 5 for the 
point of least separation av sgainst eccentricity) is given in 

figure 6. 


The fact that the shear stress shoiald decrease in the region of 
least separation of the cylinders can be seen directly from equation 
(7), which indicates that the shear stress is proportional to 

But AA.J ~ and, hence, the shear stress is approximately pro- 
portional to which decreases as the distance between cylinders 

decreases . 


The friction factor for the annulus based on the pressure gradient 
is defined by 




lb 


_ Dg dp/dx 
2P«b,av^ 

or, in terms of known dimensionless groups. 


f 




(27) 


(28) 
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PrictioG factors based on the pressiare gradient are pi^sented in figure 
7 as a function of Reynolds number and eccentricity. Comparison of these 
friction factors with those for a circular tube Indicates fair agreement 
for the case when the cylinders are concentric. The values of the fric- 
tion factor decrease as the eccentricity is increased; and, for the case 
where the cylinders touch, the friction factors are approximately 70 
percent of those for concentric cylinders. 


HEAT-TEANSFER DISTRIBUTION ON’ INNER CYLINDER 
WITH OUTER CYLINDER INSULATED 

With the information on the velocity distributions obtained in the 
preceding section and with certain assumptions, the fully developed 
heat-transfer distributions on the inner cylinder when the inner cylin- 
der is heated (or cooled) and the outer cylinder is insulated can be 
calciolated. Peripheral wall temperature distributions will be calcu- 
lated from these heat-transfer distributions in the next section. 

The heat added between two straight lines normal to the wall per 
imit length of annulus is C®®® sketch (a)') . It is assumed herein 

that all this heat is used in heating the fluid element between the 
straight lines (on both sides of the line of maximum velocities). That 
is, the net heat transfer across the sides of the element is assumed 
sm a l l compared with that transferred throu^ the wall. This assimption 
is valid when the variation of wall tenperature is not too great and the 
Prandtl number is close to or greater than unity. 

Making a heat balance on the element of fluiid between two strai^t 
lines and using the foregoing assumption result in^ 

= PS%=p(^)^ ^29) 

For the whole annulus 


p 

‘‘As an alternative assumption, an element bounded by velocity- 
gradient lines rather than by strai^t lines was used. No net heat 
transfer across the velocity-gradient lines was assumed. This assunp- 
tion was found to give essentially the same distribution of as did 

the assumption used in the text. 



12 


MGA TN 3451 


It is shown in appendix A that, for the fully developed case, 
dt-jj/dx = ^■fc'b,av/^ when the heat transfer per unit area at a given 
angle 9 does not vary with x. Division of equation (29) by equa- 
tion ( 30 ) and conversion to dimensionless form give 


(31) 

It can easily be shown that equation (31) also gives the ratio of local 
Nusselt number to average Nusselt number for the case of uniform cir- 
cumferential wall temperature, since, for that case. 


% 'C AA h 


h = 


■*^1 


Jl 

“ ’^b,av 


"*'l,av ' '*'b,av 


(32) 


and 


or 


^v 


‘^Ijav 

^l,av ~ ^,av 


( 33 ) 


h ^ Nu _ *^1 

^v ®^^v ^l,av 


( 34 ) 


Eqiiation (34) is strictly true only for the case of imiform circumfer- 
ential wall temperature, because for any other temperature distribution 
the temperature differences in equations ( 32 ) and (33) will not cancel. 

values of g.jL/ 12 . av ^^/^av uniform circumferential wall tem- 

perature) are shown as a function of angle and of eccentricity in flg- 
lire 8. As was the case with the shear stress (fig. 5), the heat trans- 
fer is low on the side of the cylinder where the separation of the cyl- 
inder is least and approaches zero when the cylinders touch. However, 
comparison of figures 5 and 8 indicates that the heat transfer ap- 
proaches zero much more rapidly than does the shear stress. !Ehe same 
trends are indicated by comparison of figures 6 and 9, where t^/t, 

and q., /q., „ ^ are plotted against eccentricity for the point of least 

separation of the cylinders. Comparison of these figures indicates, 
therefore, that the assuD 5 )tlon of similarity of the shear-stress and 
heat-transfer-coefficient distributions, which is usimlly made in anal- 
yses for odd-shaped passages, might give heat-transfer coefficients in 
the vicinity of a corner that are too high. This was pointed out in 
reference 10. Ejqperlmental wall temperature distributions in reference 
11 indicated that in some cases the predicted coefficients were too high. 
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The reason that the heat transfer decreases more rapidly than the 
shear stress as the line of least separation of the cylinders is ap- 
proached can easily he seen by coniparison of equations (?) and (29). 
Both equations contain AA/AZ]_ or which causes the heat 

transfer or shear stress to decrease near the line of least separation. 
However, the heat -transfer equation (eq. (29)) contains, in addition, 
the loc^ hulk velocity ujj, which also decreases j therefore, the heat 
transfer decreases more rapidly than the shear stress as the line of 
least separation is approached. 


WALL TEMPEMTUKE DISISIIBUTION 


The inner cylinder of the annulus is assumed, in this section, to 
he the outside surface of a thin-walled tube. Heat is transferred uni- 
formly to the inside surface of the tube. Because of tangential conduc- 
tion around the tube, the heat transfer throu^ the outer surface to the 
fliiid will not be uniform and will have the distribution obtained in the 
preceding section. Since the tube is thin-walled and there are no heat 
sources in the tube, the heat transfer per unit area throu^ the inner 
surface is equal to q^_ the average heat transfer per ■unit area 

through the ou-ter s'urface to the flioid. As in the previous section, 
the outer cylinder of the ann-ulus is insulated.. 


In order to obtain the temperature distribution around the tube, a 
heat balance is first made on an element of tube of circumferential 
length 6.1 . This heat balance gives 


^^l,av 



(35) 


where b is the thickness of the tube and q^ is the tangential heat 

conduction for unit area. Equation (35) can be -written in integral 
form as 




'^l,av^l 


“It 


(36) 


where =» 0 and q^ is zero for = 0, the point of least sepa- 

ration of the cylinders, because the -wall temperature distribution is 
symmetric about that point (the tenperature gradient is zero) . But 
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dtj_ 

“It ^ dj][ 

Substituting egiiatlon (37) In equation (36) and Integrating again resiHt 
in 



wiiere t^^^ujax wall tenrperature at the point of least separation 

of the cylinders. By using the valiies of obtained in the 

preceding section^ the difference between the maximian wall temperature 
and the wall temperature at any angle for a given heat flux can be cal- 
culated ftroni equation (38) . A dimensionless parameter containing the 
difference between the maxiTnuTn and the average wall temperatures can be 
obtained by integrating eqiiation (38), that is 




max " ^l.av' 


%,av^l 



^^^^l,max ~ ^l) 


dS 


%,avV 


(39) 


The relation for bk, (t^ - t, )/(q,T ^t^) can be obtained by sub- 

t' 1 l,av ' ' l,av 1 ' '' 

tract ing equation (38) from equation (39) . 


The results for the wall ten^jerattire distributions are presented 
in flg\are 10, irtiere bkt(t 2 _ - ■*^i,av)/C‘3i,av^l^) plotted against an- 
gle for varioTis eccentricities, and in figure 11, where 
bk^(t^^jnax ~ '^l,av^/^‘^l,av^l^^ plotted against eccentricity. It is 
of interest that t-[ - '^i^av.’ ‘'^ich is the quantity of greatest 

practical interest, can be obtained from a single curve. That is, for 
a given annulus, only two dimensionless parameters are Involved. Jor 
a given eccentricity, t-, is directly proportional to the 

uniform heat flux throtigh the inside surface of the tube. 


WaiJ. temperature distributions for varloiis other cases, such as 
the case where the tube wall is not thin or where internal heat sources 
are present in the tube wall, can also be readily calculated by using 
the heat-transfer distributions in figure 8. In order to use these 
distributions it is, of course, necessary that the outer cylinder be 
Insulated. 
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NUSSELT MJMBIE 

For obtaining all the restilts thus far (wall heat-transfer distri- 
butions, wall temperature distributions, etc.)j it was not necessary to 
use the generalized temperature distribution in figure 2. However, it 
is necessary to use that distribution for obtaining the difference be- 
tween 'VTall and b-ulk temperatures -vAiich corresponds to a given heat flow, 
that is, for obtaining the heat-transfer coefficient. 

Average heat-transfer coefficient . - The average heat-transfer co- 
efficient for the inner wall of the annulus is defined by 


l%v ~ 


^l,av 


^l,av ~ ^,av 


(40) 


where 


and 






t^u-i, dA 


av 


^,aA) 


(41) 


(42) 


(43) 


■Hie average Nusselt number corresponding to the average heat- transfer 
coefficient (eq.. (40)) can be calculated from the following equation, 
which can be verified by substituting the definitions of the various 
quantities ; 
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■where 


and 


Wuav 


^l/Pe 




% ^ 

^e 


A 


Jo 


%jav^l 


fe) 







(44) 


(45) 


(46) 


Equation (46) applies to the region between ■the inner cylinder and "the 
line of ma.xi.mum velocities. For the region between the outer cylinder 
and the line of maxiimmi velocities , y-,^ and y^ are replaced by 

and y^, respectively. As in the case of velocity distributions, the 

values of y+ are measured along normals to the inner or outer cylinder 
vhlch extend to the line of maximum velocities. The ■values of t+ are 
obtained from figure 2. In the present case, where the outer ■wall is 
Insulated, figure 2 Indicates that the temperature is uniform along the 
normals drawn from the outer cylinder to. the line of maximum velocities. 

=> 0, tg - t must also be zero, inasmuch as f*" is finite. In 

the actual case t ■woilLd ■vary some^tdiat 5J.ong the lines' normal to the 
outer cylinder. Ho^wever, a constant t is consistent ■with the assump- 
tions made in the present analysis, and the error produced is Hma.l l be- 
caxise the actual ■temperature gradients are smaJ.1 in the region between 
the outer cylinder and the line of maximum velocities. 


Average Nusselt numbers calciilated by equation (44) are plotted 
against eccentricity for a range of Reynolds numbers and for ■values 
of the parameter kr^/k^b of zero and 0.01 in figure 12. The v/all 

temperature distributions were obtained from figure 10. For the case 
where the cylinders are concentric, the Nusselt numbers are in good 
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agreement ■with "those oh'tained by the conventional Nusselt equation for 
tubes or the predicted relation for tubes in 27eference 5. The average 
Nusselt numbers decrease "with increasing eccentricity (fig. 12) as did 
the friction factors (fig. 6). For the ma-ximum eccentricity the aver- 
age Nusselt ntimbers decrease to about one-fifth -their -values for zero 
eccentricity. The curves also indicate a subs-tantlal effect of -wall 
ten 5 )erature distribution kr 2 ^/k^b on "the Nusselt numbers. Ihe Nusselt 
nxmibers decrease as kr^^/k^b Increases or as tube conductivity or 
thickness decreases. The effect of peripheral -wall temperature distri- 
bution on average Nusselt number is not generally considered in analyses 

Local heat-transfer coefficient . - The local heat-transfer coeffi- 
cient at a given angle $ is defined as 



It can be seen from equation (48) that, 
temperature, Ihi/NUg_^ » h/h^^ = ‘ll/'ll,av> 
(fig. 8). 


for uniform circumferential -wall 
as mentioned pre-viously 


Values of Nu/NUg^.^. = h/hg_^r ^or a -value of kr^^/k^-b of 0.01 are 
plotted as functions of angle 0 and of eccentricity in figure 13. 
These curves are strongly affected by changes in Reynolds number, in 
contrast to most of the preceding res-ults, because the second term in 
"the denominator of equation (48) -varies considerably -wi-bh Reynolds 
number, -whereas the first term is nearly constant. 


It ml^t be mentioned that local heat-transfer coefficients are not 
essential to give a coniplete description of -the heat transfer in eccen- 
tric annuli . Most quantities of interest co"Uld be bb-tained from the ' 
plots of -wall heat-transfer distribution, -wall -temperatxrre distribution, 
and average Nusselt nuniber. Eie local heat-transfer coefficients are 
given as a matter of interest, inasmuch as most persons concerned with 
heat transfer usually consider local heat-transfer coefficients rather 
than local heat-transfer rates. 


A 
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SUMMARY OF EISULTS 

The following resiilts were obtained from the analytical investiga- 
tion of fully developed tiorbulent heat transfer and flow in eccentric 
annuli : 


1. The maximum velocities in the annulus occurred closer to the 
inner than to the outer cylinder. The location of the line of maxi- 
mum velocities was essentially independent of Reynolds number. 

2. When the cylinders were not concentric, the velocities and shear 
stresses in the region where the separation of the cylinders was least 
were lower than the average values and approached zero when the cylin- 
ders touched. 

3. When the inner cylinder was heated and the cylinders were not 
concentric, the heat transfer to the fluid in the region idler e the 
separation of the cylinders was least was lower than the average value 
and approached zero when tlie cylinders touched. For the case idiere the 
cylinders touched, the heat transfer approached zero much more rapidly 
than did the shear stress as the line of contact was approached. 

4. The friction factors for the annulus with the cylinders concen- 
tric were very slightly higher than those for a circular tube. As the 
eccentricity increased, the friction factors decreased. 

5. When the outside surface of a thin- walled tube was assumed to 
form the inner cylinder of the annulus and heat was applied uniformly 
to the inside surface of the tube, the difference between the maximum 
and average wall temperatures for a given eccentricity was directly 
proportional to the heat flux. 

6. The average ITusselt numbers for the annulus with the cylinders 
concentric were very slightly higher than those for a tube. As the ec- 
centricity increased, the Kusselt numbers decreased. The average 
Husselt number was also found to be a ftinction of peripheral wall tem- 
perature distribution. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Cftiio, March 16,. 1955 
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APPENDIX A 


CONDITION FOR FULLY DEVELOPED EEAT TRANSFER 
For folly developed h.eat transfer and constant fluid properties^ 


■where an arbitrary heat-transfer coefficient, is independent of x. 

If hg_ were not independent of x at a great dis-fcance from, the entrance 
(cyclic variations of hg_ excluded), the absolute -value of h^ would 
become arbitrarily large as x increased so that for finite temperature 
differences the absolu-te value of q.^ would become arbitrarily large. 

The foUo-vdng equations are special cases of equation (Al) : 


■^1 “ ■feb,av 

For the case where the -wall heat transfer per unit area is independ- 
ent of X (but not of g), equations (A2) and (A3) can be differentiated 
to give the following res\ilts: 



dx dx 

dt]_ dt-[3 
dx * dx 

or 

dx dx 

The quantity d-tj^/dx is therefore Independent of 
transfer is independent of x. 


(A4) 


0 when the heat 
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APPENDIX B 


SYMBOLS 

Ihe PoUowlng symbols are used in this report: 

total area between inner cylinder and outer cyl- 
inder, sq ft 

total area between inner cylinder and line of max- 
imum velocities, sq Dt 

total area between outer cylinder and line of max- 
inmm velocities, sq ft: 

wall thickness of inner tube, ft 

specific heat of fl^Jlid at constant pressure, 
Btu/(lb)(Op) 

diameter, ft 

equivalent diameter of annulus, Dg - Dj_, ft 

conversion factor, 32.2 ft/sec^ 


local heat-transfer coefficient, 
Btu 

(sec)(sq ft)(°F) 


■^1 


Jl 

“ ^b,av^ 


average heat-transfer coefficient, 
Btu 

(sec)(sq ftj(op) 


‘^l,av 

^l,av ~ ^,av 


themal conductivity of fluid, 
Btu/(sec)(sq ft)(°P/ft) 


thermal conductivity of material of inner cylinder, 
Btu/(sec)(sq ft)(°P/ft) 


peripheral distance along wall, ft 
static pressure, Ib/sq ft abs 
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'll 


rate or heat transfer from wall of inner cylinder 
per unit area, Btu/(sec)(sq ft) 


*^l,av 


average rate of heat transfer from wall of imer 
cylinder per unit area, Btu/(sec)(sq ft) 


r radius, ft 

t temperat'ure of fluid at a point, 


arhitrairy temperature 




local hu3k temperature of fltiid at given angle at 
cross section of annulus, °F 


tQ wall temperattire, °F 

tn local wall temperature of inner cylinder or tube, 

op 


^l,max . 


maximum wall ten 5 )erature of inner cylinder or tube, 
op 


u 




X 


time-averaged velocity parallel to wall at a point, 
ft/sec 

local bulk velocity at given angle at cross section 
of annulus, ft/sec 

velocity at a point on line of maximum velocities, 
ft/sec 

axial distance along annulus, 


y normal distance from wall, ft 


^Im 

value 

of 

yi 

at 


ft 

yarn 

value 

of 

y2 

at 

•N 

II 

ft 


€ coefficient of eddy diffuslvity for momentum, 

sq ft/sec 

coefficient of eddy diffuslvity for heat, 
sq ft/sec 

9 angle made with center of inner tube, radians 
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P 

T 


absolute viscosity of fluid, (lb)(sec)/sq ft 
mass density, (lb)(sec2)/ft^ 
shear stress in fluid, Ib/sq ft 


Dimensionless quantities; 




ma x ~ ^l, av 


) 


^l,ay^l 

bkfc(ti - 

‘^l,av^l 


wall-temperature-distribution parameter where an- 
gle is zero 


wall-teraperature-distrlbution parameter 


eccentricity parameter, distance between centers 
of cylinders divided by 


friction factor. 


friction factor. 


^f^b,av 

2tt 


P%,av 


Nu 

If%v 


hDe 

ITusselt number based on local heat transfer, 

Nusselt number based on average heat transfer, 
^av^e 


n 

Pr 

Re 




constant, 0.109 
Prandtl number, c {ig/k 

ir 


Reynolds nimiber. 


pUbDg 


inner-tube radius parameter, 


4 


-r]_ dp/dx 


pyp 
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t+ 


tenrperatiire parameter. 


- ■t)cpQg'To 


U++ 


■blili-tetiqperature parameter. 


++ u-H- J A\ 


u+ 


u 


++ 


ten^eratare parameter. 


ten^ierature parameter. 


Ctj - t}CpgTj_ 


(tj_ - 


velocity parameter. 


velocity parameter. 


u 




u 




^1 


ur 


u dA 

tuik- velocity parameter, 


Un 


++ 


b,av 


balk- velocity parameter, J 


dA 




u. 


velocity parameter. 


velocity parameter. 


Ui 




u- 


Im 


vij evaluated at line of maximum velocities 
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^1 

wall' 


wall' 


yj < 


+ 

^2m 

^2 ' 


U 2 evalioated at line of maximum velocities 


a 


Subscripts i 
a 
0 
1 
2 


wall-distance parameter. 


-WF^i 

^/P ^2 


ratio of eddy diffusivity for heat transfer to 
eddy diffusivity for momentum transfer, 


/ / 

Karman constant, 0.36 


arbitrary 

pei*taining to a to, . 11 
pearbaining to inner cylinder 
pertaining to outer cylinder 
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Figure 1. - Generalized velocity distribution for fully developed turbulent flow in 
smooth tubes (ref. lO). 
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Figure 2. - Predicted generalized temperature distribution for flow of air. Prandtl number, 0.73. 
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Figure 3. - Predicted velocity distritiutions for an annulus. Reynolds number, 20,000; 
Tg^r^^, 3.5. 
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Figure 4, - 
parameter. 


(b) Reynolds number, 600,000. 

Concluded. Predlcbed variation of %/u )3 with angle 0 for various values of eccentricity 
rg/r^, 3.5. 
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Fi(Sux7 5. - Rredicted ▼«ria‘tion oi saeaar strec 

triclty parometer. 5-5^ Reynolds nun 

I 



is on inner cylajidjer irith an^e d for veurious values ox eccen- 
aber, 20,000 and. 600,000. 
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Figure 6. - Predicted variation with eccentricity parameter of 
shear stress on Inner cylinder at point of least separation. 
3.5; Reynolds numher, 20,000 and 600,000. 




Figure 7. - Predicted variation of friction factor based on static-pressure gradient vlth 

•». /-y. C 
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Fltsuie 8i - Predicted vaxiation of heat transfer from inner cylinder with angle Q for various values of eccentricity 
parameter. 3.5j Reynolds number, 20,000 and 600,000. Outer cylinder insulated. 






Eccentticity parameter, e 


Figure 9. - Predicted variation of heat transfer from inner 
cylinder at point of least separation with eccentricity 
parameter, r^/r^, 3.5; Reynolds number, 20,000 and 600,000. 

Outer cylinder insulated. 




Figure 10. - Predicted Tarlation of vail temperature on Inner tu'be vltb angle $ for rarlous ralues of 
ecceatrlcity parasketer. Constant beat flux at inside vail of inner tube; outer cylinder of annulus 
Insulated^ ^ 2 /^ 1 ^ 3.5; Reyxkolds nuznber, 20^000 and SOO^OOO. 
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Eccentricity parameter, e 


Figure 11. - Predicted variation of mEucimum teiiperature 
difference on inside tube vith eccentricity parameter. 
Constant heat flux at Inside wall of inner tube; outer 
cylinder of annulus insulated; rg/r-[_, 3.5; Reynolds 
number, 20,000 and 600,000. 
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4 



Reynolds number. Re 

^1 

(a) = 0 (uniform peripheral wall temperature dlstrihution) . 

V 

Figure 12. - Predicted variation of average Nub B elt number with Reynolds number 
and eccentricity parameter. Outer cylinder Ineulated; Prandtl number, 0.73. 
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Reynolds number. Re 

* 

=■ 0.01 (nonunifonn peripheral wall temperature distribution). 
Concluded. Predicted variation of average Nusselt number with Reynolds 


number and eccentricity parameter. Outer cylinder insulated.j Prandtl number, 0.73. 
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(a) Eeynolds nuniber, 80,000. 

Figure 13. - Predicted variation of ratio of loced hoat-tranafer coefficient to 

licftt—trazisfsr CQGfficiszili yitb 

angle 0 for varloua valuea of eccentricity parameter. Outer cylinder ineulated; r^fv-y, 3.5j kr 3 ^/k(.b, 0.01. 
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(b) Reynolds number^ 600,000. 


Figure 13 
ficlent 


. - Concluded. Predicted variation of ratio of local heat-transfer coefficient to average heat-transfer coef- 
with angle 9 for various values of eccentricity parameter. Outer cylinder insulated} rg/r,, 3.5j kr^/kj-b. 
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